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A series  of  bottom  reflection  losses  are  presented  in  graphical 
form  to  show  the  contribution  of  shear  waves  when  included  in  the 
theoretical  acoustic  model  for  bottom  sediments.  The  analysis 
illustrates  the  effect  of  introducing  shear  waves  in  the  more 
complicated  theoretical  model  for  viscoelastic  solids.  This  memo- 
randum has  been  prepared  because  it  is  believedthat  the  information 
may  be  useful  in  this  form  to  others  at  NUWC  and  to  a few  persons 
outside  NUWC.  This  memorandum  should  not  be  construed  as  a report 
since  its  only  function  is  to  present  limited  supplementary  informa- 
tion on  previously  published  theoretical  work. 
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INTRODUCTION  * 


During  the  last  twenty  years  models  of  the  bottom  sediments  for  use  in 

the  calculation  of  reflection  loss  have  become  increasingly  complex. 

The  simple  impedance  model  and  the  Rayleigh  liquid  model  frequently 

predict  inaccurate  values  of  loss.  Corrections  to  the  Rayleigh  liquid 
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model  have  been  made  by  Mackenzie,  and  Bell  and  Cole  by  assuming  the 
sound  velocity  to  be  a complex  number  with  the  imaginary  part  of  the 
velocity  related  to  sound  attenuation  in  the  sediment.  The  new  com- 
plex velocity  model  resulted  in  increased  values  of  calculated  bottom 
reflection  loss  and  better  general  agreement  with  experimental  values, 
but  neglected  the  shear  waves  that  would  be  generated  in  the  sediment. 
Since  the  attenuation  of  sound  in  sediments  was  not  accurately  known, 
the  attenuation  for  the  complex  velocity  model  could  be  chosen  for  a 
fit  of  calculated  and  experimental  values  of  "bottom  loss",  and  with- 
out regard  for  the  actual  attenuation  in  the  sediment. 

A more  accurate  model  of  the  sediments  is  that  of  a viscoelastic  mate- 
rial. Here  the  sediment  properties  are  defined  by  the  complex  Lame" 
constants.  Shear  waves  are  generated  as  a result  of  both  the  sediment 
rigidity  (related  to  the  real  part  of  p)  and  the  sediment  viscosity 
(related  to  the  imaginary  parts  of  p and  X). 

The  purpose  of  this  memorandum  is  to  compare  the  results  of  the  complex 
velocity  liquid  model  with  the  viscoelastic  model  where  the  same  density, 
compressional  velocity,  and  attenuation  in  the  sediments  are  used  for 
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both  models.  Knowledge  of  sediments  and  acoustic  measurements  sensitive 
to  the  differences  between  these  two  models  can  then  be  used  to  gain  a 
better  understanding  of  both  attenuation  and  reflection. 


* Information  from  this  memorandum  was  presented  as  Paper  T10  at  the 
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745*  Meeting  of  the  Acoustical  Society  of  America,  Miami  Beach,  Florida, 


MATHEMATICAL  ANALYSIS 
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Equations  for  the  bottom  loss  as  developed  by  Bucker  ’ are  for  a layered 

viscoelastic  solid  model.  To  compare  such  a theoretical  viscoelastic 
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model  with  work  by' Cole  and  Barnard,  et  al , the  bottom  losses  for  a 
complex  velocity  model  were  calculated  using  the  following  mathematical 
derivations  for  a water  layer  of  constant  velocity  and  an  underlying 
sediment  half-space  with  absorption. 


The  problem  of  the  reflection  of  a plane  sound  wave  at  a plane  boundary 
separating  two  media,  Fig.  A,  has  a solution  for  the  incident,  reflected, 
and  transmitted  waves  in  the  following  forms. 
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A 

exp  i 
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cos  0^  + 

k^  x 

sin  0^  - u>t) 

ref  * 
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(-  k2 
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Let  k*  = k^  sin  0^  = k2  sin  02  (by  Snell's  Law)  and  k*  a^  = k^  cos  0^ 
. . a^  - cot  . 


Then  $ inc  = A exp  i (k*  a^  z + k*  x - lot) 

$ ref  = R exp  i (-k*  a^  z + k*  x - u>t) 

$ trans  = T exp  i (k*  a2  z + k*  x - wt) . 

The  boundary  conditions  at  the  interface  are  that  there  is  continuity  of 
the  vertical  component  of  velocity  and  continuity  of  pressure. 

Thus  3^  / 3z  = 3$2  / 3z  , 

and  pj  d<P1/dt=p2  3 4>2/dt  at  z = 0 . 

Then  R = a^  p2  -a2  T = 2 a^  p^ 

al  P2  +a2  P1  * al  P2  + a2  P1 

Here  R is  the  reflection  coefficient  and  T is  the  transmission  coefficient. 

Values  of  a^  and  a2  are  determined  by  the  angle  of  the  incident  wave  and  the 
requirement  that  the  velocity  potentials  satisfy  the  wave  equation.  This 
means  that 


k*  - u sin  0^  / vQ^  , .*.  a^  ■ cot  0 . 

Also,  V $ = 1 3 , , . 

2 —2  — ,2  (the  wave  equation) 

2 3t 

Therefore,  k*^  a2  + k*^  *»  w^/v^  . 

Absorption  will  be  introduced  into  the  solution  by  making  a2  complex. 
Consider  a plane  wave  traveling  in  the  plus  z direction  in  medium  2. 
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This  wave  will  have  the  form 
exp  i (£  z - wt) 

where  £ = w/v 2 from  the  wave  equation. 

Let  ? = ?'  + i , then  the  wave  has  the  form 
exp  (i  C'  z - wt)  exp  (-  z). 

Thus  £'  is  equal  to  w/Vq2  w^ere  v^2  is  the  measured  sediment  velocity  and 
£"  is  equal  to  c^,  the  measured  attenuation  in  nepers/unit-length. 


It  follows  that 


V2  * w V02 


1 “ a2  V02 


2,22 

w + a2  vQ2 

and  a\  = c2/v2  - 1 

where  c *=  horizontal  phase  velocity  = w/k*  . 

2 

The  root  of  a2  used  in  the  calculations  corresponds  to  a wave  attenuated 
in  the  z direction. 


6 


I 


PHYSICAL  MODELS 


Three  hypothetical  physical  cases  representing  clay,  silty-sand,  and 
sand  were  considered  in  this  study.  The  same  densities,  compressional 
velocities,  and  attenuations  in  the  sediment  were  used  for  both  theoretical 
approaches,  the  complex  velocity  model  and  the  visco-elastic  model.  Bottom 
loss  curves  were  calculated  for  several  frequencies  in  the  low  kHz  region 
with  similar  results.  Only  curves  for  3.5  kHz  are  shown  in  this  memorandum. 


TABLE  Is  BOTTOM  CONSTANTS 


Velocity 

(m/sec) 

Velocity 

Ratio 

V* 

Atten(db/m) 
3.5  kHz 

Density 

S/cm 

WATER 

1530 

1.03 

CLAY 

148  A 

.970 

.84 

1.4 

SILTY- SAND 

1561 

1.02 

1.8 

1.7 

SAND 

1701 

1.12 

1.116 

1.9 

sediment  velocity 


SHEAR  WAVES 


Shear  waves  have  been  found  to  exist  in  marine  sediments.  Shear  waves 
were  measured  "in  situ"  by  Bucker  in  shallow  water  and  to  3500  ft.  depth 
on  the  continental  shelf.  He  found  measurable  shear  velocities  in  mud 
as  well  as  sand.  High  shear  velocities  of  approximately  260  meters/sec 
for  deep  sea  clays  have  been  measured  by  Hamilton  ^ on  stations  of 
two  recent  cruises  in  the  Pacific.  To  correspond  to  a rigorous  theore- 
tical examination  of  bottom  reflection  loss,  shear  waves  would  necessarily 
need  to  be  included. 


Rigidity  of  the  sediment  is  introduced  in  the  theoretical  viscoelastic 
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solid  model  by  the  ratio,  Re/i/ReA  . This  is  equivalent  to  specifying 
Poisson's  ratio  for  the  sediment  layer.  It  should  be  noted  that  when 
Re^i/ReX  =0,  some  effect  of  the  shear  waves  is  still  present  from  the 
contribution  of  Im^t  related  to  viscosity.  Cases  were  computed  for 
ratios  of  0,  0.01,  0.03,  0.1,  and  0.2.  These  ratios  are  equivalent 
to  specifying  a range  of  Poisson's  ratio  from  0.5  to  0.417. 


• RESULTS 

The  bottom  loss  was  calculated  for  two  models,  both  of  which  consisted 
of  a water  layer  overlying  a sediment  half-space.  Velocity  and  attenua- 


tion selections  varied  the  sediment  from  clay  to  silty-sand  to  sand. 
Varying  amounts  of  rigidity,  r = Re^/ReX  , were  introduced  into  the 
sediment  layer.  The  results  are  shown  in  figures  1 to  3.  As  seen  in 


The  dashed  lines  labeled  "no  shear"  represent  the  complex  velocity, 
or  absorbing  liquid,  model.  When  the  ratio  r is  zero  there  is  still 
discrepancy  between  the  two  curves,  indicating  the  contribution  of  the 


imaginary  part  of  n , which  may  be  significant  even  in  this  case 
when  the  two  models  are  most  nearly  alike.  As  the  number  of  sediment 
layers  increases  beyond  the  single  sediment  half-space  a more  decisive 
difference  than  illustrated  here  would  be  expected  in  the  results. 

Figure  1 shows  bottom  reflection  loss  plotted  versus  grazing  angle 
for  a low-velocity  bottom.  Both  models  result  in  "intromission"  losses 
and  there  is  a shifting  of  the  high-loss  angular  region  as  "r"  is  in- 
creased in  the  visco-elastic  model.  Note,  however,  that  the  result 
for  r = 0.03  would  be  similar  to  the  "no  shear  wave"  model.  Plots 
for  r = 0.1  and  0.2  are  not  realistic  for  clay  but  are  included  here 
for  completeness. 


Figure  2 is  a silty-sand  case.  The  sound  velocity  ratio  of  sound 

velocities  in  sediment  and  water  for  n slightly  > 1,  the  attenuation 

and  density  are  typical  of  silty-sand  bottom  sediments.  Increasing  the 

r ratio  greatly  increased  the  bottom  loss  from  20  to  30  degrees.  From 
7 

sea-floor  data  , r ratios  from  0.1  to  0.2  seem  to  be  reasonable  for 
sediments  of  this  velocity  and  density. 

Figure  3 represents  a "fast  bottom"  case  of  sand.  Notice  a change  of 
scale  on  bottom  loss.  With  no  attenuation  one  would  have  the  Rayleigh 
case.  The  Rayleigh  form  seems  recognizable  here,  but  with  an  added  loss 
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as  the  critical  angle  is  approached.  This  probably  indicates  that 
attenuation  is  not  as  important  here  to  the  shape  of  the  curve  as  for 
the  low-velocity  case  in  Figure  1. 

The  "r"  ratios  from  measurements  tend  to  run  0.2  and  above.  In  this 
case,  r = 0.03  and  0.1  should  not  be  considered  realistic;  and  even 
higher  values  of  bottom  loss  might  be  calculated  and  observed.  The 
difference  in  models  near  40  and  50  degrees  seems  to  be  measurable  and 
suggests  further  experiments. 

Figures  4 and  5 illustrate  changes  of  attenuation  a in  sediments. 

Both  treat  a constant  "low  velocity"  bottom  where  rj  is  slightly  < | 
which  is  appropriate  for  some  clays. 

Figure  4 shows  the  effect  on  bottom  loss  for  the  complex  velocity  model 
when  the  velocity  was  held  constant  and  the  attenuation  values  were 
varied  from  0 to  10  db/meter.  There  was  less  bottom  loss  with  increased 
attenuation.  It  should  be  noted  that  for  low-velocity  liquid  sediments 
the  angle  of  intromission  occurs  when  a p,  = a o as  can  be  seen 

from  the  equation  for  the  reflection  coefficient.  For  sediments  with 
attenuation  this  condition  cannot  be  satisfied  as  will  be  complex. 

Figure  5 shows  the  effect  of  the  same  increases  of  attenuation  for  the 
visco-elastic  model  when  r = 0.  As  the  sediment  attenuation  increases, 
the  bottom  losses  at  first  become  less  near  the  angle  of  intromission  In 
this  model.  However,  with  further  increasing  attenuation  in  this  visco- 
elastic model  a peak  in  the  bottom  loss  curve  will  eventually  build  up 
as  shown  for  a = 10.  This  reappearance  of  a peak  loss  is  due  to  the 
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attenuation  of  shear  waves  in  the  sediment.  Although  the  real  part  of 
ti  = 0,  there  is  a contribution  from  the  imaginary  part  of  n . This 
peak,  though  occurring  at  a very  high  attenuation  for  our  purpose,  does 
show  the  effect  of  shear  attenuation  in  a highly  viscous  material, 
i 

CONCLUSION 

* 

An  effort  was  made  to  select  three  representative  sea-floor  types. 

In  all  cases  the  introduction  of  shear  waves  resulted  in  higher- theore- 
tical bottom  reflection  losses.  For  the  case  of  the  shear-wave  model 
with  zero  rigidity  bottom  reflection  loss  was  very  close  to  the  "no 
shear"  model.  For  non-zero  rigidity  and  reasonable  ranges  of  "r" 
ratios  we  feel  that  the  results  definitely  are  not  the  same  for  cases 
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of  interest  and  for  the  two  types  of  models.  It  is  believed  that  at 
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least  a few  important  sea-floor  areas  do  exhibit  these  rigidity  ratios, 
and  that  for  a rigorous,  accurate  treatment  the  shear  waves  should  be 
included  in  the  theoretical  model. 

Also,  these  calculations  imply  that  for  normal  incidence  reflection 
there  is  not  much  difference  for  either  model  or  for  any  of  the  sedi- 
ments investigated.  For  sands,  and  grazing  angles  30  to  50  degrees, 
there  should  be  measurable  losses  that  might  be  used  to  demonstrate 
differences  between  models. 
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